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We consider the propagation of strong gravitational waves interacting with a nonperturbative vac¬ 
uum of spinor fields. To described the latter, we suggest an approximate model. The corresponding 
Einstein equation has the form of the Schrodinger equation. Its gravitational-wave solution is anal¬ 
ogous to the solution of the Schrodinger equation for an electron moving in a periodic potential. 
The general solution for the periodic gravitational waves is found. The analog of the Kronig-Penney 
model for gravitational waves is considered. It is shown that the suggested gravitational-wave model 
permits the existence of weak electric charge and current densities concomitant with the gravitational 
wave. Based on this observation, a possible experimental verihcation of the model is suggested. 
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I. INTRODUCTION 

An experimental search for gravitational waves (GW) is one of the most intriguing problems in modern physics. 
The discovery of GWs will give us confidence that we are moving in the right direction in understanding classical 
gravity. If we cannot detect GWs, then either our apparatus is not sensitive enough or something is not taken into 
account when considering the propagation of GWs. 

One of the effects accompanying the propagation of GWs could be its interaction with a nonperturbative spinor 
vacuum. The reason for such an interaction to occur is that the energy-momentum tensor of a spinor field contains 
the spin connection, which in turn contains first derivatives of tetrad components with respect to the coordinates. 
As a result, the Einstein equations give the wave equation for a GW which contains second derivatives of the tetrad 
components on the left-hand side and their first derivatives on the right-hand side. 

In Ref. [I| we have considered the propagation of a weak GW interacting with the nonperturbative spinor vacuum. 
Here we extend those results to the case of a strong GW. In doing so, as in Ref. [I|, to model the nonperturbative 
vacuum of a spinor field, we will use a phenomenological approach. Within the framework of this approach, we 
make some physically reasonable assumptions about expectation values of the spinor field and its dispersion. This 
will permit us to reduce the infinite system of differential equations for all Green’s functions of the nonperturbative 
quantum spinor field to the finite set of equations (for more details, see Refs. Qi). 

Following this approach, here we will discuss the solution of the Einstein equations for a strong GW propagating 
on the background of the nonperturbative vacuum of spinor fields, which is a generalization of the weak, plane 
gravitational wave of Ref. [I| . 


II. STRONG GW IN A NONPERTURBATIVE SPINOR VACUUM 


According to the textbook Q, the metric for a strong GW propagating 


in one direction is sought in the form 


ds^ = 2d^dri gab{ri)dx°‘dx^, 


( 1 ) 
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where =77 are lightlike coordinates, and the indices a, b run over the values 2, 3. It is convenient to introduce 

the new variables 


9ab{r]) = -X^{i)lab{ri), (2) 

det(7ob) = 1. (3) 

The only nonvanishing component of the Einstein tensor G^u is (/I is the tetrad index and v is the spacetime 
index), so that one can obtain the following Einstein equation for the strong GW propagating in a flat vacuum 
spacetime: 


X + ^ X = 0. (4) 

Here the dot denotes differentiation with respect to rj; 7 “^ is the two-dimensional tensor reciprocal to jab- The function 
xiv) is an unknown function, and jab are arbitrary functions, obeying the constraint ([3]). 

As pointed out in Ref. the presence of the term jacj^‘^jbdj°'^ has the result that after a finite time interval 

X becomes zero. This in turn leads to vanishing of the metric determinant g, i.e., a singularity in the metric. But 

this singularity is not physically significant since it is related only to the unsatisfactory nature of the reference frame, 
“spoiled” by the passing gravitational wave, and can be eliminated by a suitable coordinate transformation. 

Our goal here is to consider the propagation of GWs in a nonperturbative spinor vacuum. We expect that such a 
physical system has to be considered in nonperturbative language when both a metric and a spinor field are regarded 
as quantum quantities and are quantized in a nonperturbative manner. Then, to describe the interaction between a 
quantum metric and a quantum spinor field, we have to write down the Einstein-Dirac operator equations. 

Our nonperturbative approach for quantizing nonlinear fields is described in Ref. [l[ . Within this method, we have, 
strictly speaking, to solve the following operator equations: 

Gatj. = xfatj., (5) 

^'4’ — m'4 = 0 , ( 6 ) 

where G^a is the operator of the Einstein tensor; T^i, is the operator of the energy-momentum tensor; ip is the operator 
of the spinor field; a is the tetrad index; 77 is the coordinate index; is the covariant derivative for the spinor with 
the appropriate spin connection; = SttG, G is the gravitational constant. Hereafter we use units where c = h = 1. 

As mentioned in Ref. [l[, the solution for such a set of operator equations can be found by writing down and 
subsequent solving an infinite system of equations for all Green’s functions for quantum fields. In practice, such a 
procedure can be carried out only approximately. This means that we have to use some physical arguments to cut off 
the above-mentioned infinite set of equations, and then solve the resulting finite set of equations. 

Our approach is based on the following assumptions: 

(i) The metric remains always classical. 

(ii) The spinor field is quantum. 

(iii) The spinor field is decomposed as a product of q— and c—numbers. 

(iv) In order to calculate the energy-momentum tensor of the spinor field, we assume some ansatze for a 2-point 
Green’s function of ip- 

In such an approximation the Einstein-Dirac operator equations (O and (| 6 ]) can be written in the following man¬ 
ner [l|: 


Q 


\7 T-^ 


Gafi — ^ \ Q 

q) = 0 


T-a 


Q), 


( 7 ) 

( 8 ) 


where (<5 | - | <5) is the quantum averaging with respect to the quantum state | Q). 

Also, to check the consistency of the ansatze for a 2-point Green’s function, instead of solving the Dirac equation (jS]), 
we will use the Bianchi identities (|S]), as we did in Ref. [l|. 
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III. APPROXIMATE MODEL OF THE NONPERTURBATIVE SPINOR VACUUM 


In order to write down the energy-momentum tensor for a vacuum of a spinor field interacting with a GW, we must 
have a model of a nonperturbative vacuum of the spinor field. Let us emphasize once more that we cannot use a 
perturbative model of a spinor field since in the presence of a gravitational field the set of equations m and m is a 
strongly nonlinear system. 

Taking this into account and following Ref. [I| , our strategy in the formulation of the model of the nonperturbative 
spinor vacuum is as follows: (i) we write some classical ansatze for a spinor; (ii) we derive the corresponding energy- 
momentum tensor; and (iii) we then write hats over the corresponding spinor components. 

As the first step, we take the ansatze for the classical spinor in the form 


i) = 



(9) 


For the tetrad p3l) below, there is only one nonvanishing component Consequently, we have to choose the 

components of the spinor A, B, U, Q so that we obtain the corresponding component Ti^. This happens \i B = A and 
Q = V. In this case we have 


Ti^ = —2VV* (/3' cosh (3 + a' sinh /3 -|- 4a;), 


( 10 ) 


where the prime denotes d/drj. 

With the ansatze Q and the energy-momentum tensor (HOD in hand, we suggest the following approximate model 
of the nonperturbative vacuum of the spinor field. 

• The nonperturbative vacuum is described by the following operator of the spinor field: 




A 

V 

\vJ 


( 11 ) 


The constant operators A, V appearing here are independent of rj. 
• The corresponding energy-momentum tensor of the spinor field is 


Q 


Ti 


Iri 


q'^ = —2 iP' cosh/3 -I- a' sinh/3 -I- 4a;). 


( 12 ) 


To check this model, we calculate the divergence of the energy-momentum tensor and show that it vanishes. 


IV. EINSTEIN EQUATIONS FOR A STRONG GW INTERACTING WITH THE NONPERTURBATIVE 

SPINOR VACUUM 


For the metric o, we seek a solution of the Einstein equations for a strong GW propagating through the nonper¬ 
turbative spinor vacuum. The tetrad for this metric is 


/O 1 0 0 \ 

10 0 0 

0 0 x(7/)e“(’')/2 0 

Vo 0 sinh/3(?7) 


(13) 


with the corresponding two-dimensional metric 


^ab — 


Here a{r]) and /3(?7) are arbitrary functions Q. 


e“ cosh^ /3 sinh /3 
sinh /3 


(14) 
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The substitution of the tetrad (fT3l) and the energy-momentum tensor (fT^ into the Einstein equations ([71) yields 
the equation 


■ X 


/?' cosh (i + a' sinh /3 ~ ^ ) cosh^ ^ 


X = -4wx, 


(15) 


where the prime denotes differentiation with respect to the dimensionless fj = k (vv^j r]^ and the dimensionless 

ib = Loj (For convenience, we omit the tilde from fj and u) hereafter.) The terms /?'cosh/3, a'sinh/3 on 

the left-hand side and 4a; on the right-hand side of Eq. (11511 are the imprints of the nonperturbative spinor vacuum. 
One sees immediately that Eq. (fTKll is a Schrddinger-like equation with the effective potential 

14//=/3'cosh/3-I-a'sinh,5 — i -b/3' ^ cosh^/3 — sinh2/3. (16) 

In what follows we will seek periodic solutions to Eq. dTK]) . It is clear that for periodic functions a,/3 we obtain an 
equation similar to that describing the movement of a single electron in a crystal. Such an equation has been well 
studied in the literature (see, for example, the textbook [^), and we can apply all mathematical methods used in 
solving the Schrodinger equation to our case of a strong GW. 


V. GENERAL SOLUTION 

In order to solve Eq. (fTSl) for a periodic potential, in this section we apply the methods developed in solid state 
theory. The only difference is that the function xiv) is a real function, unlike the usual quantum mechanics where a 
wave function is complex. 

First of all, let us rewrite Eq. (US in the form 

- x"(»7) + yeff{v)xiv) = (17) 

In order to find a solution of the type of GW, we have to investigate the case of a periodic potential 

yeff{v + Vo) = Veff{v), (18) 

where rjo is the period of the effective potential. 

We seek a solution to Eq. (flTl) in the form 

OO 

xiv) =Xo + Y^ 

k=l 

a set of equations for the coefficients xo,cik,bk, we: (a) integrate Eq. (ITTl) : (b) multiply (flTl) by 
integrate; and (c) multiply dig by ^ sin and integrate. This yields the following system 

of equations: 


In order to obtain 
— cos and 

Vo V Vo ') 


( 2'Kk 

ak cos - r] 

\ VO 


( 2'Kk 

bk sin - r] 

\ VO 


(19) 


-I OO 

14+2 ^2 4 ” bkVk 

^ k=l 

2 00 

) XO^n ^ ^ ^k^nk 

00 

bn + Xoin + [^afcl4ri + 


27rn 

Vo 


k=l 


= -4a;xo, 


= —4a;a„ 


= —4a;6„ 


( 20 ) 

( 21 ) 

( 22 ) 
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where n = 1,2, ■ ■ ■ ,oo and 


+>?o/2 


Vh = — / Veffdv; 

Vo J 


Vk = 


Vk = 


^nk — 


^nk — 


^nk — 


-i)o/2 

+170/2 

2 f /27rfc \ „ 

— / 14//cos - V ] dv, k = 1,2,-■ ■ , 00 ; 

VO J V ^0 y 

-rial‘s 

+170/2 

2 , (2'Kk \ , , , „ 

— / Ve//sin - ?7 077, fc = 1,2,--- , 00 ; 

VO i \ VO ) 

-170/2 

+170/2 

2 f f 27rn \ / 27rfc \ 

— / cos -?7 Vef f COS - rj drj, n, k = 1,2, ■ ■ ■ , 00 ; 

Vo J \ Vo J \ Vo J 

-vo/2 


2 

Vo 


+170/2 


f / 2TTn \ / 27rfc \ 

/ cos -77 14// sin -77 077, n, fc = 1,2, • • • , 00; 

J \ Vo J \ Vo J 

-vo/2 
+170/2 

2 f ■ f 27r7T, \ , / 27rfc \ 

— / sin -77 14// sin -77 (277 n, fe = 1, 2, • • • , 00. 

Vo J \ Vo J \ Vo J 


(23) 

(24) 

(25) 

(26) 

(27) 

(28) 


In principle, by using the set of equations one can find a GW solution for any periodic metric functions 

a( 77 ) and P{v)- Remarkably, in doing so, one can employ all well-developed methods of solid state theory (see, for 
example, the textbook i). 


VI. THE KRONIG-PENNEY MODEL FOR GRAVITATIONAL WAVES 


In the previous section we have used the well-known methods of solid state theory to obtain the general solution 
for the GW. To discuss some particular features of the GW, here we employ an approximate approach of solid state 
theory (the Kronig-Penney model) for the description of the GW. 

Let us consider the simplest case, /3 = 0, for which Eq. dni yields 


// 1 /2 . 

- X - X = -dwy. 

For such a case we have the following two-dimensional metric (fHl) : 


Vab — 


0 


0 e- 


(29) 


(30) 


In order to apply the Kronig-Penney approximation, we will use the metric function 0 ( 77 ) in the range —a < rj < b: 

(31) 


Go 

0 ( 77 ) = <( Vb?? + Gi 


if — a < T] < 0, 
b 

2 ’ 


if 0 < ?7 < |, 

b 
2 


-Vov + C 2 if I < 77 < 6 , 


where Go, 1,2 and Vq are constants; a -I- I* = 770 , and 

a{r] + nr]o) =a{r]), 77 = ±1, ±2, • • • . 

The corresponding effective potential 14// = — is 


^eff = 


0 if — a < 77 < 0 , 
—if 0 < r] < b 


(32) 


(33) 
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with the periodicity 


Veff{r] + nr]Q) = Veff{r]), n = ±1, ±2, • • • . (34) 

The profiles of a and 14 // are shown in Fig. [TJ 

To solve Eq. with the effective potential ( 1551 ) , we will employ the Kronig-Penney method exactly in the same 
way as in solid state theory but remembering that we have to seek a real solution instead of a complex one, as is done 
for the Schrodinger equation. 

For the periodic potential ( 1551 ) . Eq. (|2^ takes the form 


I x" = K'^x if - a < 7? < 0, 

\x" = -Q^X if 0<V<b, 

where = 4a;, = 14^/4 — K^. 

For the fundamental region —a<rj<b, we seek a solution in the following form 

+ if -a<77<0, 

1 c sin Qrj + d cos Qp if 0 < rj < b. 

Using these, we first constrain continuity inside the fundamental domain (i.e., for jy = 0): 


(35) 


(36) 


\f + 9 = d, 

[Kif -g) = Qc. 

For the region b < rj < a + b, we seek a solution in the form 

X{p) = 

We then constrain continuity outside the fundamental domain (i.e., for ry = 6): 

j csinQb + d cos Qb = fe~^°‘ + ge^‘^ , 

1 cQ cos Qb — dQsinQb = fKe~^°‘ — gKe^°‘. 


(37) 


(38) 


(39) 


A necessary condition for the set of equations ( 1571 ) and (I39|) to have a nontrivial solution is that the determinant of 
the corresponding matrix is zero: 


0 

Q 

sin Qb 
Q cos Qb 


1 -1 -1 

0 -K K 

cos Qb — 

-QsinQ6 -7X6"^“ 


(40) 


The resulting constraint equation for the parameters K,Q,a, and b is 




Q^ 


sinh aK sin bQ + cosh aK cos bQ = 1. 


2KQ 


(41) 
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FIG. 2: The profile of the left-hand side of the constraint equation (143 II for a = 3 and b = 50. The straight line corresponds 
to 1. 


Introducing the dimensionless parameters 


X = 




b 



(42) 


we obtain from ((411) 


2a;- 1 
2Vx — x' 


: sinh (a^yiPj sin (i\/l — x^ + cosh (di/x) cos (bVT 


= 1 . 


(43) 


Here 0 < a; < 1, i.e., 0 < w < V)f/16. The latter expression is obtained from the condition > 0 which has been 
assumed to be satisfied in the derivation of the above results. The case of < 0 is discussed below. 

The typical profile of the left-hand side (Ihs) of the constraint equation (H51) is shown in Fig. [2] One can see that 
there is a finite number n of possible values of which depends on the values of the parameters a and h. The 
numerical values of Un are determined by the points of intersection of the lines in Fig. [21 Note that in the limit a; —>■ 0 
the Ihs = —(d/2) sin5 -(- cos 5, and when a; = 1, the Ihs = (&/2) sinhd -I- coshd, i.e., the Ihs remains always finite. 

The case of < 0 can be obtained from (|41l) by changing Q to iQ'. 


K^ + Q^ 
2KQ 


sinh aK sinh bQ -|- cosh aK cosh bQ = 1. 


The corresponding dimensionless equation will then be 
2a;- 1 


2y/x‘^ — : 


sinh (di/s) sinh (i^/x^^^ + cosh (di/x) cosh (b^/3^ 


- 1 = 1 . 


(44) 


(45) 


Here a; > 1, i.e., to > V)//16. One can easily see that the Ihs of this expression is always greater than 1, i.e., Eq. (|45ll 
does not have a solution. Physically, this means that there are no GWs with w > Vj//16. 


VII. EXPERIMENTAL VERIFICATION 

Let us discuss now possible experimental consequences coming from a consideration of the propagation of a strong 
GW through the nonperturbative spinor vacuum. 















Direct calculations show that for such a case there exists the following electric current 


= 


(4(0 


VV^ 


g),0,0,0 


( 46 ) 


where = ^,r],y,z. 

Making the transformation to Cartesian coordinates t, x as 


^ t — X t -\- X 

f ^ = 






and similarly, 


we have 




At _ AX 

ji = 1 _ I- jV = _ 

^ V2 ' ^ 72 ’ 


f = -f = 2V2(^Q\vV^\Qy 


(47) 


(48) 


(49) 


This means that we have the electric charge and current densities concomitant with the GW. The electric current is 
directed along the direction of propagation of the GW. This observation allows us to suggest the following experimental 
verification of the GW and nonperturbative spinor vacuum models considered here: It is possible to try to measure 
the weak electric charge and current densities together with the standard measurements of GWs (LIGO, LISA, and so 
on). On the other hand, instead of measurements of the electric current, one can measure the corresponding magnetic 
field. 


VIII. THE CORRECTNESS OF THE ANSATZE FOR THE SPINOR FIELD 

The correctness of Eq. dH) is verified directly. To do this, one can calculate the expression V as the classical 
one and then substitute in it the 2-point Green’s function Then, taking into account that we have only one 

nonvanishing component of the energy-momentum tenor [see Eq. pOI) ]. the Bianchi identities ([8]) are trivially 
satisfied. 


IX. CONCLUSION 


To summarize our results: 

• We have suggested an approximate model of the nonperturbative spinor vacuum. 

• The propagation of strong gravitational waves interacting with such a vacuum has been investigated. 

• It was shown that the Einstein equation reduces to a Schrodinger-like equation with a periodic potential. 

• We have obtained the solution for the special case of the metric function corresponding to a GW with the gyy,zz 
components only. 

• A possible experimental verification of the nonperturbative spinor vacuum model interacting with a strong 
gravitational wave has been suggested. 

Eor better qualitative understanding of some features of the resulting gravitational waves, we have considered the 
simplest case where the periodic potential is modeled by the Kronig-Penney potential of solid state theory. As a 
result, it was shown that the GW parameters (the wave length 770 and the amplitude Vq) can be arbitrary but the 
parameter affecting the form of the plane wave in the spinor vacuum, is quantized. 

The correctness of the approximate model of the nonperturbative spinor vacuum was verihed by calculating the 
divergence of the right-hand side of the Einstein equations. One might expect that for an improved model of the 
nonperturbative vacuum it will be necessary to use the nonperturbatively quantized Dirac equation where the term 
“spin connection x spinor field” (i.e., {Q \u}abfj,'4t\ Q)) need to be taken into account. The latter means that we have 
to take into account the quantum correlation between a metric and a spinor field. 
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Lastly, it may be noted that the formal analogy of the GW equation here obtained with the Schrddinger equation 
allows the possibility of using the well-developed methods from other fields of physics to solve such an equation. In 
particular, for the periodic metric functions discussed here, Eq. (Ull) is identical in form to the equation of motion 
of an electron in a crystal. Aside from the formal aspect, such an analogy would lead us to use physical intuition to 
obtain further results in this area, as it is done, for example, in analogue gravity 


Acknowledgments 

This work was partially supported by a grant in fundamental research in natural sciences by the Science Committee 
of the Ministry of Education and Science of Kazakhstan and by a grant of the Volkswagen Foundation. 


[1] V. Dzhunushaliev and V. Folomeev, Eur. Pliys. J. C 74, no. 9, 3057 (2014) arXiv:1405.6121 [gr-qc]]. 

[2] V. Dzhunushaliev, V. Folomeev, B. Kleihaus and J. Kunz, Eur. Phys. J. C 74, 2743 (2014) arXiv:1312.0225 [gr-qc]]. 

[3] L. Landau and E. Lifshitz, The Classical Theory of Fields (Pergamon, Oxford, 1987). 

[4] W. A. Harrison, Solid State Theory (Dover Publications, 2011). 

[5] C. Barcelo, S. Liberati and M. Visser, Living Rev. Rel. 8, 12 (2005) [Living Rev. Rel. 14, 3 (2011)] gr-qc/0505065 . 



